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Semi-Annual Status Report on 
The Structure of Evaporating and Combusting Sprays: 
Measurements and Predictions 


SUMMARY 

This report covers the first semi-annual period of an investigation 
of the properties of noncombusting and combusting sprays. Activities 
during this period were limited to noncombusting sprays. An apparatus 
was constructed to provide measurements in open sprays with no zones of 
recirculation, in order to provide well-defined conditions for use in 
evaluating spray models. The measurements in noncombusting sprays in- 
clude: mean and fluctuating gas velocities, mean and fluctuating drop 
velocities, drop sizes, mean liquid fluxes, mean temperatures and mean 
compositions. Thus far, measurements have been completed in a gas Jet, 
in order to test experimental methods, and are currently in progress for 
nonevaporating sprays. 

Models of the process are also being examined including: (1) a locally 
homogeneous flow (LHF) model where interphase transport rates are assumed 
to be infinitely fast; (2) a separated flow (SF) model which allows for 
finite interphase transport races but neglects effects of turbulent fluctu- 
ations on drop motion; and (3) a stochastic SF model which considers effects 
of turbulent fluctuations on drop motion. In this report, the models are 
evaluated using existing data on particle-laden Jets. The LHF model gener- 
ally overestimated rates of particle dispersion while the SF model under- 
estimated dispersion rates. In contrast, the stochastic SF flow yielded 
satisfactory predictions except at high particle mass loadings where effects 
of turbulence modulation may have caused the model to overestimate turbu- 
lence levels — and thus dispersion rates. This evaluation is tentative, 
however, due to uncertainties in the specification of initial conditions for 
the existing data base. 

Current work is addressing both existing and present new measurements 
in nonevaporating sprays. These results should indicate the importance of 
drop collisions and shattering, as well as effects of turbulence modulation. 
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Introduction 


Th« potential value of rational design procedures for llquld~fuelad 
coabustora has motivated continuous efforts to develop reliable aodele of 
spray evaporation and combustion processes. The goal la to reduce the 
time and cost of cut and try methods of development by providing a better 
understanding of fundamental spray proceasSs and methods for estimating 
the effect of specific design changes. While there ha^ been significant 
progress in the development of spray models, no existing model has demon- 
strated a capability for reliable predictions (or even correlation) of 
spray behavior, laro.»ly due to the unavailability of appropriate measure- 
ments for model evaluation. The objective of this investigation is to 
complete measurements in evaporating and combusting sprays in order to 
close this gap in the existing information. 

Current activity concerning spray modeling is quite high, as evidenced 
by Che recent development of relatively sophisticated models of sprays in 
industry [l“5l. Reviews of spray modeling [6-?l] find 10-20 papers per 
year reporting new spray models at present, when contributions from 
government and universities are considered. While the basic elements of 
spray models are in hand, much work must still be done to improve the treat- 
ment of interphase transport, turbulence, chemical reaction and radiation 
processes in sprays [6-11]. 

Properly establishing models by comparing predictions with experi- 
ments has substantially lagged model development. The major problem la 
that measurements in sprays, particularly combusting sprays, present 
formidable experimental difficulties. A second problem is chat experi- 
menters have not adequately recognized the need for specifying the proper- 
ties of the injector. Spray models are so dependent on Injected drop sizes 
and velocities that measurements withe.. ^ this information are of little 
value for testing spray models. 

The objective of this investigation is to complete new measurements 
in order to improve the data base available for model development. A main 
feature of the experiments is tiuit a steady axisymmetric spray, with no 
zones of recirculation. Is considered. Such conditions minimize experi- 
mental difficulties, reduce uncertainties in modeling the turbulent flow, 
and simplify requirements for obtaining accurate numerical solutions of 
Che transport equatiens. Analysis of the experimental results is also 
being undertaken using models whose development was initiated during past 
NASA-sponsored research in this laboratory. 

The investigation is planned for a two-year period, with noncombusting 
sprays examined during the first year and combusting sprays during the 
second year. This report covers the first semi-annual period of the in- 
vestigation where attention was primarily devoted to analysis of particle- 
laden flows and measurements in nonevaporating sprays. This report 
considers experimental and theoretical methods used during the investi- 
gation, along with a description of the results obtained to date. A paper 
describing a portion of these results was also prepared during this report 
period, cf. Ref. 13. 


Numbers in brackets denote references. 
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2< Experlaantal Methods 

2.1 T rat Apparatus 

A akatch of tha apray apparatus appaars In Fig. 1. Tha praaant 
nonconbusting flows all hava dansiciaa graatar than air; tharaforsi tha 
Injactor is diractad downward in still air» in ordar to avoid racirculation. 
Tha maasuranants anploy optical instruaantation which sust ba Mountad on a 
rigid basa. Tharafora, probing the flow is acconplishad by travarsing 
tha injactor in threa dimanaions. 

The flow is protactad froa room disturbancas using a scraanad 
enclosura a squara by 2.5 a high). Kajor travarsing, to obtain radial 
profiles of flow quantities, involves aoving tha antira caga asaaably. 

This keep^i the flow nearly concentric with the vertical axis of tha caga, 
whicii siniaizes disturbance of tha axisyaaetric flow due o off-cantar 
screen positions. 

The inlet to the exhaust system is screened and is located 1 a below 
the plane of instrumentation. Testing has shown that operation of tha 
exiuiust system has a negligible effect on flow properties at tha measuring 
plane . 

A Spraying Systems Co. air atomizing injector (model 1/4J2050 fluid 
nozzle and 671A7 air nozzle with outlet diameter of 1.19 mm) is being 
used for the spray tests — similar to earlier work {12,13]. An attempt 
was made to employ a flashing injector developed earlier [16]; however, 
this arrangement proved to be unsatisfactory due to complications of 
external expansion from a choked injector for conditions yielding desired 
atomization properties. 

The flow system of the injector is illustrated in Fig. 2. Filtered 
dry air is supplied from a storage tank. The flow rate of air is con- 
trolled with a pressure regulator and metered with a critical-flow orifice 
to ensure long-term stable operation. 

The liquid is stored in a tank under pressure, however, the tank is 
not agitated and pressure levels are moderate (0.3-0. 8 kPa) . Therefore, 
the dissolved air content of the liquid is negligible. The flow rate of 
liquid is controlled with a valve and metered with a rotameter. In order 
to maintain repeatable flow and atomization conditions the entire test 
cell is heated above normal ambient temperatures - to 27 ^1 C. 

The test liquid for measurements with nonevaporating sprays is vacuum 
pump oil. This material has reasonable pumping cliaracteristics and very 
low volatility. Freon 11 will be used for tests of evaporating noncom- 
busting sprays, since It has sufficient volatility to evaporate rapidly 
in room temperature air while posing no explosion hazards. Freon 11 was 
also used in previous experiments with sprays in this laboratory [12]. 

2.2 Instrumentation 


Table 1 is a summary of ttie measurements that are being undertaken 
in noncombusting sprays and the techniques employed for each measurement. 
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Kic. 1 Sketch of tl»e experimental apparatus. 



CRITICAL FLOW 



Sketch of the injector flow system 








Table 1. Summary of Instrumentation for Noncombusting Sprays 
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A laser Doppler anemometer is used to measure mean and fluctuating 
gas velocities. Several beam orientations provide measuremeiits of 
various velocity components as well as the Reynolds stress. In order 
to avoid problems of fringe bias and reverse flows» the laser beams are 
frequency shifted. Concentration biasing and effects of drops are 
avoided by employing high concentrations of seeding particles— following 
past methods [12»13]. 

Accurate measurement of initial conditions requires good spatial 
resolution in the dense portion of the spray. Imaging techniques are 
most effective in this region [17]. Therefore, drop properties are 
determined by high-magnification, direct-lighted streak photographs 
near the injector. This method is tedious, however, it is the most 
r'*liable technique available for providing initial conditions rela- 
tively near the injector. 

A second method used for drop size and velocity measurements Involves 
the LDA — visibility technique [18]. This technique is relatively inde- 
pendent of the refractive indices of the disperse and continuous phases 
and can be used in noncombusting sprays. The technique provides a 
simultaneous indication of drop size and velocity far more conveniently 
than the imaging method — in dilute sprays. This equipment is cur- 
rently being acquired — funded by University resources — and the system 
is not yet operational. However, we hope to complete a portion of the 
noncombusting spray measurements during the next report period with this 
instrument. 

Drop sizes are also measured with a Malvern particle size rnalyzer, 
which operates on the principle of Fraunhofer diffraction of laser light 
scattered by drops. This instrument only provides a linc-of-sight mea- 
surement and is largely used to monitor injector performance. 

Liquid fluxes are measured by isokinetic sampling and by impaction. 
These techniques are well-establisho.d in the literature and provide a 
check of the newer optical methods. The sampling probe has a 2 mm ID 
sampling port, followed by a nitrogen-flushed diverging section to prevent 
drop impaction on the internal surfaces of the probe. The drops are col- 
lected In an internal spill-over probe, backed up by a filter to collect 
small-diameter drops. Isokinetic sampling conditions are determined by 
fixing the probe just below the LDA measuring volume and matching dis- 
turbed and undisturbed gas velocities. 

It is difficult to obtain accurate mean temperature measurements in 
sprays due to effects of drop impaction. We have developed a shielded 
temperature probe which circumvents most of these problems fl2,13]. This 
involves a 50 ym diameter Pt/Ft 10% Rh butt-welded thermocoaple located 
along the centerline of a 1.6 mm OD semi-circular shield. The shield is 
positioned upstream of the probe so that all but the smallest drops 
(diameter less than '0 ym generally) impinge on it. Therefore, the 
thermocouple primar. senses gas temperatures. The probe has been 
checked in gas flames and found to yield temperatures within 10 K of 
unshielded temperatures, indicating that effects of flow deflection and 
increas'^d heat losses are small for current experimental conditions. 
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The mean composition of Injected fluid in the evaporating apraya ia 
determined by Isokinetic sampling, at the mean gas velocity, and analyais 
with a gas chromatograph. The probe is heated to prevent blockage %rlth 
liquid. This tec ique provides a direct measurement of the local mixture 
fraction of the flow. The system used is identical to Shearer et al. [12]. 

3. Theory 

3.1 General Description 

Present measurements as well as results in the literature are 
being used to evaluate spray models. Two models are considered: (1) a 
locally homogeneous flow (LHF) model, where both phases are assumed to 
have the same velocity and^temperature and are in thermodynamic equilibrium 
at each point in the flow; and (Z) a separated flow (SF) model where finite 
interphase transport rates are considered, so that each phase generally has 
a different velocity and temperature throughout the flow. 

Both models employ the widely adopted procedures of k-c-g turbulence 
models for the gas phase [12-15,19,20], since this approach has been 
thoroughly calibrated during earlier work in this laboratory [12-15]. 

Major assumptions for the gas phase are: exchange coefficients of all 
species and heat are the same, buoyancy only affects the mean flow, and 
kinetic energy and heat losses are negligible. Effects of buoyancy and 
radiation heat losses are generally small in practical sprays; therefore, 
treating these phenomena as perturbations is Justified. Neglecting 
kinetic energy limits the model to low Mach number flows, which is 
appropriate for the test conditions to be examined, as well as most 
practical combustion chambers. The assumption of equal exchange coef- 
ficients is widely recognized as being appropriate for high Reynolds number 
turbulent flows typical of spray processes. 

Two versions of the separated flow model have been developed. One 
version considers mean particle trajectories similar to past models while 
the second version treats effects of turbulent fluctuations using stochastic 
methods. 

In order to ensure adequate numerical closure with reasonable compu- 
tation costs, the model is limited to boundary-layer flows with no recircu- 
lation. The present test flows are axlsymmetric with no swirl; therefore, 
the analysis is posed accordingly. An advantage of these conditions is 
that they correspond to cases where k-c-g models were developed and have 
high reliability. 

3 . 2 Locally Homogeneous Flow Model 

The governing equations, under the LHF approximation, can be 
written as follows: 


* 

This approximation is only valid for infinite Interphase transport 
rates — implying a flow with Infinitely small particles or drops. 
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The parameters (|) and S. appearing in Eqs. (1) and (2) are summarized in 
Table 2, along with the appropriate empirical constants. The empirical 
constants were established for noncombusting single-phase flows and are 
not changed when sprays are considered [12,13>20]. The turbulent vis- 
cosity is calculated from the turbulent kinetic energy and rate of dissi- 
pation as follows: 


Pt • C^Pk^/e (3) 

With equal exchange coefficients and moderate radiative heat losses, 
which is normally the case even for combusting sprays. Instantaneous 
properties are only a function of mixture fraction (essentially the 
state resulting from adiabatic equilibration of f kg of Injector fluid 
and (l-f)fg of ambient flu-^d) . This allows determination of all prop- 
erties, temperature, composition and density, as a function of mixture 
fraction — once and for all. Properties in the flow are then found from 
the probability density function for mixture fraction, P(f) as follows 
[ 20 ] 




f 4>(f)P(f) df 


(A) 


where 4>(f) is known from the state relationship. A functional form must 
be assumed for P(f). A clipped Gaussian function, similar to Lockwood 
and Naguib [20], is being used for P(f), although the specific form used 
has little effect on predictions. P(f) depends on two parameters which 
can be found from the local values of f and g as described by Lockwood and 
Naguib [20]. 

The computation of the state relationships, (j)(f) involves 
only conventional adiabatic mixing (noncombusting spray) 
or adiabatic flame (combusting spray) calculations. 

The methods used for these computations are described elsewhere 
[11-14,21,22]. In the case of combusting sprays, the state relationship 
calculations employ CEC-72 [23], which allows for variable specific heats 
and dissociation in high-temperature gases. 

A typical state relationship is illustrated in Fig. 3. The case 
shown is for a Freon 11 spray injected into air at atmospheric pressure, 
using an air-atomizing injector [12]. Evaporation of Freon 11 in the 
atomizing air yields a mixture of air and Freon-11 liquid and vapor 
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Table 2. Source Terns 
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leaving the injector at a temperature somewhat below their initial 300 K 
temperature (at f «• 1) . The surroundings are pure air at 300 K (at f ■ 0) . 
As fluid mixes with the injector flow, the concentration of liquid and the 
temperature decrease until the point is reached where all the liquid has 
evaporated (f ~ 0.3) » which also corresponds to the maximum Freon 11 vapor 
concentration. At lower mixture fractions, the temperature increases and 
the Freon 11 vapor concentration decreases toward the ambient values. 

The presence of liquid in portions of the mixture fraction rang'' causes a 
very nonlinear density variation. Similar state relationships ur com- 
busting flows at moderate and high pressures are presented elsewhere 
[11,13,141. 

Initial conditions for the calculations are prescribed at the injector 
exit similar to past work [11-14]. In the absence of other information, 
all properties are assumed to be constant at the injector exit, aside from 
a shear layer having a thickness equal to 1 percent of the injector radius 
at the passage wall. The constant property portion of the flow is 
specified as follows: 


X “ 0, r < 0.99d/2; u “ M /m , f 

o o o o 


1, k 

o 


(0.02G 

o 


g 


o 


0 


e = 2.84 X 10‘^ u ^/d 
o o 


(5) 


The Initial variation of properties in the shear layer is assumed to be 
linear. Equation (5) provides the inner boundary condition until the 
shear layer reaches the spray axis, after which all gradients at the axis 
are zero. The ambient values of u, f, k, c and g are all zero for the 
flows to be considered here. The calculations are performed using a 
modified version of GENMIX [24]. 

The main advantage of the LHF model is that there are only a few 
empirical constants, which are specified by earlier measurements, and 
only routine equilibrium computations and simplified injector quantities 
are required for computations. In view of its ease of use, the model 
has demonstrated reasonably good capabilities for computing spray 
properties [12-14]. The main defect is that the rate of flow develop- 
ment is overestimated in the two-phase region. Potential users of spray 
models are likely to begin with a version of this type, therefore, it is 
desirable to evaluate its performance. 

3.3 Separated Flow Model 

The separated flow model adopts the main features of the LHF 
model, but only for the gas phase. The liquid phase is treated by 
solving the Lagrangian equations of motion of the drops and then computing 
the source terms in the governing equations for the gas phase, resulting 
from interphase transport processes. This general procedure corresponds 
to the particle tracking or particle-source-in-cell methods used in most 
recent two-phase models and is often called a discrete droplet model (DDM) 
[ 11 ]. 
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Drop Model . The main assumptions of the drop trajectory calcu- 
lations are as follows: dilute spray with drop transport parameters 
equivalent to a single drop in an Infinite environment} ambient con- 
ditions given by mean flow properties; negligible effect of turbulent 
fluctuations on drop transport rates; empirical treatment of drag and 
convection effects; quaslsteady gas phase; negligible drop shattering 
and collisions; liquid surface in thermodynamic equilibrium; and 
negligible radiation, Dufour and Soret effects. These assumptions are 
conanon in spray models— their Justification is discussed elsewhere (7,11]. 

The assumption of unity Lewis number is not made, since liquid fuels 
are often modeled poorly in this manner due to their large molecular 
weights. Average gas-phase properties for drop transport calculations 
are found at a reference state defined by 


4 * 


avg 


■a 4> + (1 - a) i 

%x> 


( 6 ) 


where is a generic property representing temperature and species mass 
fractions. The factor u is selected to match measured transport rates 
during supported drop calibration experiments. The value a ■ 0.9 has 
generally been found to provide the best fit of the data during past work 
[11,12], Mixing rules and values used for property calculations are 
described elsewliere [11,12], 

Due to the difficulty of completely modeling Internal transport 
processes of drops, the analysis considers two limits: (1) the infinite 
liquid diffusivity approximation where all properties within a drop are 
assumed to be uniform at each instant of time; and (2) the negligible 
diffusivity or "onion skin" model, where the drop surface adapts im- 
mediately to changes in local ambient conditions, while the bulk liquid 
remains at its initial state. These cases bound the range exhibited by 
real drops in sprays. 

Drops leaving the injector are divided into n groups, defined by the 
initial position, sise, velocity and direction. The subsequent motion of 
each group is given by 


pi 


X + 

poi 


u .dt 

pi 


(7) 


The instantaneous velocity is determined by integrating the equation for 
conservation of momentum 


du , 
dt 


- ("/«) it 


^ , -V 


pi 


- ‘ u) + a 


( 8 ) 
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where the site and mass of the drop are related as follows: 

- (TT/6) jdJ^ (9) 

assualng uniform liquid properties. Conservation of nass of the drop 
yields 



( 10 ) 


where m is determined from the local concentrations at the liquid surface 
and in the gas [11,12]. 


Conservation of energy at the drop, for tne infinite liquid dlffu~ 
sivity model, yields 




dT , 




( 11 ) 


for a simple one component liquid. In the case of the onion skin model 
there is no bulk heating; therefore, the energy equation becomes 


TTd ^ h (f-T ) - m (h - h ,) (12) 

p ' s ' sg pf 

where T is not equal to the bulk liquid temperature T * T any longer. 

In both^cases, the gas phase mass fraction of drop liqSid is given by the 
vapor pressure relationship. The present calculations are limited to low 
pressures, the Clausius-Clapeyron equation suffices and the concentration 
of dissolved gas is negligible. The general vapor pressure relationship is 





(13) 


The coacentratlons of other species at the drop surface, needed for 
property evaluations, are found from film theory analysis [7,11] 


isg i* Fsg 


(19) 


Drop transport is modeled using film theory. This involves analysis of 
transport rates at Re <■ 0 and then multiplying the resulting expression by 
a correction for forced convection [7,U]. The basic expression for the 
drop mass transfer rate is 
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\e-0 

27rd PD 
P 


In ( 




) 


while the expression for the heat transfer coefficient is 


\e-od Ac Ac 


^ird X 
P 


'2TTd X' 
P 


(16) 


The multiplicative correction for convection was established by Faeth 
and Lazar [25] 


h A 

or — 


I + 0.278 Re^ ^ ^ CPr or Sc)^ ^^ 


\e-0 ®Re-0 (1 + 1.232/ (Re (Pr or 


(17) 


The standard drag coefficient for solid spheres is employed in the calcu*> 
lations, approximated as follows [15]: 


2/3 


(18) 


0.44 


, Re > 1000 


Gas-Phase Model . This portion of the analysis utilizes the 
dilute spray approximation. This Implies that the void fraction is unity 
and that the presence of drops does not contribute directly to the gener- 
ation or dissipation of turbulence (the latter effect is frequently called 
turbulence modulation) . 


The interaction between the liquid and gas phases yields an additional 
source term on the RUS of Eq. (1). These terms are found by computing 
the net change of mass and momentum of each drop class i as it passes 
across a computational cell J. The source terms are as follows: 


Mass: 


Momentum; 


n 


• “ 1 n, (m. . “ m. ^ ^) , 

dmj i i in i out j 


(19) 


'iZi ■ Ji *l“"l“pl>out - <Vpl'ln>j 


( 20 ) 


Only one index is needed to define a cell since the calculation Is 
parabolic and each computational cell is defined by its radial position. 
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page is 
quality 

Mixture Fraction: n^(n^ (2X) 

where A Is the number of drops per unit time In each class, in and out 
denote ^rop conditions entering and leaving the computational cell, and 
Y . denotes the mass of liquid per unit mass flow rate of injected fluid. 
tHIs procedure allows for full Interaction between tlM phases, which is 
vital for treating the near injector region. 

The gas-phase equations are solved in the same manner as the LHF 
model. The only change in this portion of the program involves addition 
of the new source terms given by Eqs. (13)-(21). The particle motion 
equations, Eqs. (7)-(18), are solved at the same time, in a stepwise 
fashion, using a second-order finite difference algorithm. The solution 
requires the specification of the properties of all drop classes at fhe 
Injector exit. 

The properties of the gas phase are computed from f and g in the same 
manner as the LHF model. The main change in this portion of the calculation 
is that the state relationships, and f, only pertain to vaporized liquid. 
This procedure Is simplified by Ignoring energy absorbed by drop heating 
in the infinite liquid dlffuslvlty model. This effect could be handled 
by solving an additional equation for conservation of energy. in the gas 
phase; however, the error is small at moderate pressures and the added 
complexity of the formulation is not merited at present. 

3.4 Stochastic Separated Flow Model 

The basic SF analysis provides only for deterministic trajectories 
for particle groups. In practical turbulent flows, however, particles are 
also dispersed by turbulent fluctuations. The SF model was extended to 
treat this effect, since the phenomena can be an Important mechanism for 
the spread of drops in sprays. The approach used to describe particle dis- 
persion adapts a stochastic technique proposed by Gosman and loannldes [26]. 

The stochastic model Involves computing the trajectories of a statis- 
tically significant sample of Individual particles as they move away from 
the injector and encounter a random distribution of turbulent eddies — 
utilizing Monte Carlo methods. The general nature of the stochastic 
approach is illustrated in Fig. 4 for a single particle encountering three 
eddies. The key elements of this approach are the methods for specifying 
the properties of each eddy and the time of interaction of a particle with 
a particular eddy. The k-e-g representation of turbulence used in the SF 
model provides a convenient method for prescribing these quantities. 

Properties within a particular eddy are assumed to be uniform, but 
properties change in a random fashion from eddy to eddy. The trajectory 
calculation is the same as the basic SF model, involving solution of 
Eqs. (7)-(18); however, mean gas properties in these equations are re- 
placed by the instantaneous properties of each eddy. 
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Sketch of the stochastic particle trajectory model. 
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The properties of each eddy are found at the start of Interaction by 
making a random selection from the pdf of velocity and mixture fraction. 

The velocity fluctuations are assumed to be isotropic with a Gaussian ... 
probability density distribution having a standard deviation of (2k/3) ' . 
The cumulative distribution function for the three velocity components 
Is formed and each distribution Is randomly sampled. This Involves 
selecting three numbers In the range 0-1.0 and computing the velocity 
components at these three values of the cumulative distribution function. 

Instantaneous physical properties for each eddy are found In a similar 
manner. The Instantaneous mixture fraction ls_assumed to have clipped 
Gaussian probability function with mean value £ and variance g — which 
completely specify the pdf of mixture fraction, cf. Lockwood and Naguib 
[20]. The cumulative distribution function Is constructed for this pdf 
and a single random number selection in the range 0-1.0 yields the 
Instantaneous mixture fraction of the eddy at the sample value of the 
cumulative distribution function. The state relationships then provide 
the physical properties of the eddy, at this mixture fraction, as 
described earlier. 

A particle Is assumed to Interact with an eddy for a time which Is 
the minimum of either the eddy lifetime or the transit time required for 
the particle to cross the eddy. These times are estimated assuming that 
the characteristic size of an eddy Is the dissipation length scale. 


L 

e 




( 22 ) 


The eddy lifetime Is computed from the following: 

t - L /(2k/3)^^^ 
e e 


(23) 


The transit time of a particle is found from the linearized equation of 
particle motion in a uniform flow 

t- - X In (1 - L /(Tiu'’-u"|)) (24) 

t e ' p 

where u" - u” is the velocity difference at the start of the interaction 
and X Is the**linearized particle relaxation time 


X 


dp/(3pc„|u" 



(25) 


When L > x|u" - u"|, the linearized stopping distance of the particle is 
smallef than the cRaracterlstlc length scale of the eddy and Eq. (24) has 
no solution. In this case, the eddy has captured the particle and the 
Interaction time Is the eddy lifetime. 


*These prescriptions are slightly different than those of Gosman and 

loannldes [26], but have been successfully calibrated [15]. 
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The remainder of the computation proceeds similar to the deterministic 
case. The only change Is that the source terms of Eqs. (19) **(21) are com- 
puted for the random walk trajectories of the particles— as opposed to 
their deterministic solution. The main disadvantage of the stochastic 
method Is that more particle trajectories must be considered In order to 
obtain statistically significant particle properties. 

The stochastic model yields estimates of both mean and fluctuating 
particle properties at each point In the flow. This information is useful » 
since these properties can be measured and provide a good test of model 
performance. A notable feature of the model is that added empiricism Is 
minimal . ** 

4. Results and Discussion 


The apparatus was assembled during this report period and measure- 
ments have begun. The model was evaluated using existing results for 
particle-laden flows during the report period. 

In the following, results for particle-laden flows will be considered 
first, followed by the Initial findings for nonevaporating sprays. 

4.1 Particle-Laden Flows 


A search of the literature revealed a substantial data base of 
measurements In particle-laden Jets. Table 3 is a summary of the results 
employed to evaluate the models. The theoretical findings of Hlnze [27], 
concerning the dispersion of infinitely small particles from a point 
source In an Isotropic turbulent flow were employed to calibrate the 
stochastic model, similar to Gosroan and loannides [26]. The remaining 
measurements were employed for model evaluation. This Included measure- 
ments of single-particle dispersion in a turbulent duct flow by Snyder 
and Lumley [28]; and in particle-laden jets by McComb and Salih [30], 

Laats and Frishman [31], and Levy and Lockwood [32]. 

The calibration results are illustrated in Fig. 5. Stochastic mode;, 
predictions from Gosman and loannides [26] and the present model are shown, 
along with the exact analytical results from Hlnze [27]. The present model 
is in good agreement with the exact results, indicating satisfactory cali- 
bration. Gosman and loannides [26] model, as posed, yields poorer agree- 
ment due to a computational error in their original work that has since 
been corrected [33]. 

Predictions of the present stochastic SF model, as well as the version 
of Gosman and loannides [26] and the measurements of Snyder and Lumley 
[28] are Illustrated in Fig. 6. These experiments involved dispersion of 
individual particles which were isoklnetlcally injected into a uniform 
turbulent flow downstream of a grid. The two models yield nearly the 


In fact, no new constants must be formally prescribed; however, Eqs. 
(22) and (23) implicitly involve constants of proportionality although their 
values turned out to be unity when the model was calibrated. 


Table 3. Sunnary of Measurements In Particle-Laden Flows 
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Fig. 5 Analytical and stochastic solutions for the 
dispersion of small particles in homogeneous 
isotropic turbulent flow with long diffusion timts. 
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Fig. 6. Predicted and measured particle dispersion in a uniform, grid-generated turbulent flow. 
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same results for heavy particles where the calibration error of [26] has 
little effect. The present model yields better agreement for the lighter 
particles where calibration of characteristic eddy time and length scales 
is more significant since light particles tend to be captured by the 
eddies. 

Significance of particle dispersion in particle-laden Jets can be 
better appreciated from the results appearing in Fig. 7. The particle 
concentration measurements of Yuu et al. [29] are Illustrated along with 
three predictions: (1) the LHF model, (2) the SF model with no consider- 
ation of turbulent dispersion, and (3) the stochastic SF model. Pre- 
dictions in each case are given at values of x/d which limit the data 
range. The LHF model overestimates the rate of particle spread, since 
it neglects slip between the phases. On the other hand, the deterministic 
SF model underestimates the rate of spread, since radial particle vel- 
ocities are only generated by the v velocity component — which is rela- 
tively small (this version also Indicates that the particles are confined 
to a progressively narrower range of r/x as x/d increases, which is also 
not observed). In contrast, the stochastic SF model yields good agreement 
with the measurements. 

Predictions of the LHF and stochastic SF models are compared with 
data from Yuu et al. [29] at various operating conditions in Fig. 8. The 
stochastic SF model yields good predictions while the LHF model continues 
to overestimate the rate of particle spread. 

The measurements of McComb and Salih [30] are compared with the model 
in Figs. 9 and 10. The results are generally similar to those obtained 
with the Yuu et al. [29] measurements. An interesting feature of these 
measurements is that the particle spread rate begins to approach the LHF 
predictions as x/d becomes large, since the large eddies in this region 
tend to capture the particles more effectively. 

Predictions and measurements for the test conditions of Laats and 
Frishman [31] are Illustrated in Figs. 11 and 12. These authors considered 
much higher mass leadings of particles (mass of particles/mass of gas) than 
Refs. 27-30. The agreement between the stochastic SF predictions and mea- 
surements is comparable to that observed for the Yuu et al. [29] and McComb 
and Salih [30] data, except at the highest particle mass loading where 
the stochastic model overestimates the rate of spread of the particles. 

This suggests that turbulence levels are lower under these conditions 
than predicted, which could result from the extraction of turbulence 
energy by the drag of the particles. This phenomena, which is called 
turbulence modulation [34], is not considered in the present models. 

Additional evidence of effects of particles on turbulence properties 
Is provided by the results of Levy and Lockwood [32] illustrated in 
Fig. 13. Mean gas velocities are predicted reasonably well; however, 
longitudinal fluctuations of gas velocity are generally underestimated 
for large particles at high loading ratios. In this case, the high rate 
of slip of the particle through the gas enhances turbulent fluctuations 
in a manner that Is not considered in the model. 
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Fig, 7 Turbulent particle dispersion in an 

axisynunetric jet Data of Yuu et al. [29], 
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Predicted and measured axial gas velocities and particle 
concentrations in a dust-laden air Jet with particle 
diameter of 20 ym. Data of Yuu et al. [291. 
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Fig. 10 Predicted and measured axial velocities and particle concentrations 
in particle-laden jets. Data of McComb and Salih [30]. 
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FIr* Predicted and measured axial mean velocities along the centerline 
of particle-laden Jets. Data of Laats and Frishman (31]. 
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Fig. 13 Predicted and measured mean and fluctuating gas velocities 
in particle-laden jets. Data of Levy and Lockwood (32]. 
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Figure 14 is an illustration of mean and fluctuating particle 
velocities for the data of Levy and Lockwood [32]. Predictions of mean 
particle velocity agree reasonably well with measurements at all con- 
ditions. Predicted particle fluctuating velocities, however, under- 
estimate the measurements for all but the smallest particles in this 
data base. This effect can only be attributed to cperatlon of the 
particle-gas mixing system (a worm gear drive followed by a straight 
but relatively short length of ducting) since the large particles have 
a relatively moderate change in velocity over the flow field. Unfor- 
tunately, Levy and Lockwood [32] did not make measurements of initial 
conditions at the duct exit to help clarify this effect. 

Our analysis of existing data for particle-laden Jet was inhibited 
by lack of information concerning Initial conditions in all cases. The 
original authors generally felt that there was negligible slip between 
the phases at the Injector exit; however, analysis of the injector 
passage flow invariably showed that this was not the case. Thus, «re 
conclude that the results of the model are encouraging; however, more 
Information on initial conditions is needed for decisive eval\iatlon. 
Furthermore, the model may require extension to include effects of 
turbulence modulation before adequate predictions can be obtained in 
dense particle flows and sprays. 

4.2 Nonevaporating Sprays 

The apparatus described in Section 2 was designed, fabricated 
and assembled. Testing thus far has been devoted to the isothermal air 
jet and noncombusting sprays. 

The main objective of the air jet experiments is to check the oper- 
ation of the apparatus and to establish LDA procedures for measuring mean 
and fluctuating gas velocities. The new results are being compared with 
existing measurements of Shearer et al. [21], Mao et al. [22], Wygnanski 
and Fiedler [36], and Uetsronl and Sokolov [37]. The measurements involve 
axial profiles of mean velocity and streamwise velocity fluctuations as 
well as radial profiles of these quantities and Reynolds stress. These 
measurements are being conducted for a subsonic alr^jet passing from the 
injector, with jet Reynolds numbers greater than 10 , to ensure reasonably 
well-developed turbulent flow. Jet exit conditions are defined by mea- 
suring the air mass flow rate and the thrust on the Injector. 

The measurements are being conducted with the single-channel LDA, 
using various beam orientations to obtain the velocity components and 
the Reynolds stress — similar to earlier work [21-22]. Results obtained 
to date indicate that adequate seeding levels have been achieved, flow 
disturbances are negligible, the jet is properly aligned, and the results 
of the measurements are in good agreement with past work. 

Testing has begun to examine the properties of nonevaporating sprays. 
Thus far, two conditions have been established for detailed structure 
measurements: (1) a spray with a SMD of 30 microns which represents a 

near LHF condition, and (2) as spray with a SMD of 70 microns which re- 
quires SF modeling for reasonable accuracy. Thus far initial conditions 
have been measured — SMD, mass flow rates and injector thrust. Detailed 
structure measurements are now in progress using the LDA for gas velocities. 
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Modeling efforts ere now being directed no noneveporetlng sprays. 

LHF aodel calculations have been completed for present test conditions. 

SF model calculations will be undertaken as soon as Initial conditions 
for the experiments have been fully defined. In addition to these mea- 
surements. recent results reported by Yule et al. [35] are also being 
considered. These results with nonevaporating sprays will help indicate 
whether modeling of steady spray processes will require consideration of 
particle collisions and shattering. 

5. Status and Plans for the Next Report Period 

The spray apparatus has been constructed and measurements are cur- 
rently In progress for nonevaporating sprays. Based on progress to date, 
all testing on noncombusting sprays should be complete at the end of the 
current grant period. The extent to which the LOA visibility method for 
measuring drop-size and velocity distributions can be applied to these 
results will depend on prompt delivery of equipment and problems encoun- 
tered in making the system operational. The results of the Investigation 
are not contingent on the use of this technique, however, since the 
availability of the equipment was not known at the time the study was 
planned . 

The LHF and SF models have been developed and coded for use on the 
computer. Turbulence model constants used for these calculations were 
established during earlier work [21-22]. Computations for solid particle 
flows have generally established the usefulness of the stochastic particle 
dispersion model for dilute flows, although some uncertainty remains due 
to Inadequate specification of initial conditions throughout the existing 
data base. Results to date suggest chat there may be significant effects 
of turbulence modulation In dense particle flows, however, additional In- 
formation on flow structure Is needed to rationally complete this extension. 
Comparisons between predictions and both exis.:ing [35] and the present new 
measurements will be undertaken during the next report period. These re- 
sults will help determine whether model extensions are needed In order to 
treat drop collisions and shattering, as well as effects of turbulence 
modulation. 
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